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We investigate the effect of spatially correlated disorder on two-dimensional topological insulators 
and on the quantum spin Hall effect which the helical edge states in these systems give rise to. 
Our work expands the scope of previous investigations which found that uncorrelated disorder can 
induce a non-trivial phase called the topological Anderson insulator (TAI). In extension of these 
studies, we find that spatial correlations in the disorder can entirely suppress the emergence of the 
TAI phase. We explain this phenomenon by generalizing an existing effective medium theory to the 
case of correlated disorder potentials. The predictions of this theory are in good agreement with 
our numerics and may be crucial for future experiments. 
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I. INTRODUCTION 



A two-dimensional topological insulatoi^ features edge 
states similar to those of the quantum Hall effect with 
the difference that electrons of different spin move in the 
same direction at opposite edges. Accordingly, this so- 
called "quantum spin Hall effect" (QSHE)2 can be un- 
derstood as two non-interacting copies of a quantum Hall 
system, one for each spin. The stability of the edge states 
is guaranteed by time-reversal symmetry which forbids 
scattering into the counter-propagating edge state with 
opposite spin.^El These properties, which have recently 
attracted considerable attentionpHin] ma k; e topological in- 
sulators promising cand idate s for key components in fu- 
ture spin-tronic devices™^. i n 2 006 HgTe/CdTe quan- 
tum wells were proposed as suitable systems for a first 
experimental realization of the QSHE 13 which was, in- 
deed, achieved shortly thereafter JSGH Numerically it was 
found that the edge states in a topological insulator not 
only show great robustness against disorder but also 
that strong disorder itself can induce a phase featur- 
ing pure edge transport even if the system is an ordi- 
nary insulator in the clean limit This disorder-induced 
and topologically non-trivial phase was named topolog- 
ical Anderson insulator (TAI). In 2009 a theory was 
put forward^ that lead to a detailed understanding of 
the TAI, showing that disorder causes a negative cor- 
rection to the topological mass which pushes the sys- 
tem into the TAI phase. This interesting phenomenon 
has meanwhile been investigated numerically in a va- 
riety of different systems^HZl] including the case of a 
three-dimensional topological insulator! 2 ^ However, due 
to the challenges involved in controlling the disorder in a 
HgTe/CdTe quantum well, the TAI has not yet been re- 
alized experimentally. This problem might be overcome 
by employing ultra-cold atomic gases in optical lattices 
for the realization of a topological insulatorP^l In such a 
highly tunable model system the disorder could be intro- 
duced by an optical laser speckle potentiaP^H which has 
the advantage of being under external control. 
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FIG. 1. (color online). Illustration of the scattering setup for 
the considered topological insulator: a rectangular disordered 
middle part of height L y and width L x is attached to two semi- 
infinite leads on the left and right. Random potentials with 
and without correlations between neighboring grid points in 
the underlying lattice discretization are shown in the bottom 
and top panel, respectively. 



An important point to emphasize in this context is that 
both the speckle pattern for cold atomic gases as well as 
the disorder which naturally occurs in a quantum well 
are characterized by a finite spatial correlation length £. 
Since this correlation has been disregarded in all pre- 
v ious numeric al studies of the TAI which we are aware 
ji6 | i7 | 22 | 27 | 28 | ^ ne q ues ti n was posecPSl how a finite cor- 
relation length £ would influence the predictions for the 
appearance and for the stability of the TAI. In view of 
the fact that spatial correlations in the disorder have al- 
ready been shown to play an important role in the con- 
text of various other scattering scenarios^HUl one may 
expect such correlations to be a relevant factor also for 
topological insulators. We address this topic by studying 
explicitly how a static and spatially correlated disorder 
influences the transport characteristics of topological in- 
sulators (Fig [I]). As we will specify in detail below, our 
numerical results show marked deviations from conven- 
tional simulations with uncorrelated disorder. 
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II. METHOD 



T at the Fermi energy Ep. 



We proceed along the lines of previous studies, 
where an appropriate description of the two dimensional 
HgTe/CdTe quantum well was proposed in terms of an 
effective HamiltonianP^ This Hamiltonian which was de- 
rived based on the k • p perturbation theory and the six- 
band Kane-model takes the following form: 
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and a 1 labelling the Pauli-matrices. The basis of this 
effective Hamiltonian consists of the s-like El and the 
p-like heavy-hole HI quantum well sub-bands for spin 
up (+) and down (— ). The ordering is chosen to 
be \El+),\Hl+), \E1-), \H1-). Since the spin-up and 
spin-down parts in the Hamiltonian are decoupled as 
a consequence of time reversal symmetry^ it is suffi- 
cient for our calculations to only use the spin-up block 
h(k). The solution for the spin-down block follows from 
a time-reversal operation. The material-dependent con- 
stants A, B and D in all our calculations are set to re- 
alistic values A = 364.5 meVnm, B — —686 meVnm 2 , 
D = —512 meV nm 2 and C = meV taken from Ref. 
[44] . The sign of the topological mass m has a strong 
impact on the system's transport behavior: for positive 
in the system behaves like an ordinary insulator with a 
band gap of 2 |m|, whereas if m is set to a negative value 
the system turns into a topological insulator featuring 
perfectly transmitting edge states for the Fermi energy 
Ep lying inside the bulk band gap \Ep\ < |m|. To simu- 
late such a system we use the experimentally determined 
valued of m — — 10 meV. 

The scattering geometry which we consider consists of 
a rectangular disordered region of width L x and height L y 
attached to two clean, semi-infite and undoped leads (see 
illustration in Fig. [I]). We discretize the scattering region 
on a square lattice with grid constant a using the tight- 
binding approximation in the continuum limit for the im- 
plementation of the effective Hamiltonian. If not stated 
otherwise, the grid-constant a is set to 5 n m, in agree- 
ment with the value used in previous studies F^ 6 * 17 * 27 * 45 *^ 
For simplicity, we consider the limit of vanishing tem- 
perature T — > and infinitely small bias voltages V — » 
applied between the two semi-infinite leads. According to 
the Landauer-Buttiker formalism the conductance G in 
this limit is proportional to the transmission probability 
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The indices n and m extend over all N lead modes 
and t nm labels the transmission amplitude from mode 
n in the incoming lead to mode m in the outgoing 
lead. Since we consider both spins separately, every 
mode only contributes a single conductance quantum 
e 2 /h. For the calculation of the transmission we em- 
ploy the advanced modular Recursive Green's function 
methocP^^ w hich incorporates the disorder by way of a 
static on-site energy-value V(x) imposed at every grid- 
point x = (xi,yj). In previous calculations the random 
on-site energies V(x) were chosen to be uniformly dis- 
tributed within a given energy interval [— U/2,U/2] and 
each random sample from this distribution was drawn in- 
dependently for each grid point. Since in this case the 
disorder value on each grid point has no correlation with 
the values on neighboring grid points we will refer to 
this type of disorder as "uncorrelated" . To go beyond 
this limitation and to account for the spatial correlations 
which naturally occur in realistic situations we choose 
our disorder potential such as to obey the Gaussian cor- 
relation function 



C(r) = (V(x) ■ V(x + r)) cx exp - 
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where the brackets (...) stand for an average over all the 
grid points x and many disorder realizations. The stan- 
dard deviation of this Gaussian defines the correlation 
length £, which measures the spatial range of the correla- 
tions. The value of the disorder strength U is established 
by demanding 
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These values are chosen such as to agree with those of 
the uncorrelated disorder potential distributed within the 
interval [-U/2, U/2]. See Fig. [ljfor an illustration of the 
disorder potentials with and without spatial correlations 
in the employed tight-binding grid. 



III. RESULTS 

We first consider the conductance through a rectangu- 
lar bar of a topological insulator (width L x = 2000 nm 
and height L y = 500 nm) with topological mass m ~ 
— 10 meV. In the clean limit (see U = in Fig. [2]ja)) 
quantized edge transport (green area) only occurs in an 
energy range \Ep\ < |m|. 16 Adding now an uncorrelated 
disorder to the clean topological insulator extends this 
QSH phase and shifts its onset towards higher energies 
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(see U > in Fig. ga) and Ref. [16]). This so-called 
TAI phase of quantized transport emerges in the presence 
of strong uncorrelated disorder and extends beyond the 
limits \Ep\ = \m\. Our results from Fig. [2ja) agree very 
well with the literature}^ thereby confirming the validity 
of our simulations. 

In a next step we repeat this calculation for a correlated 
disorder potential. We choose the value of the correla- 
tion length £ = 23.45 nm considerably larger than the 
grid constant a = 5 nm but still much smaller than the 
height L y = 500 nm of the sample. Our results for such 
a finite correlation length (see Fig. |2^b)) differ dramati- 
cally from the uncorrelated case (see Fig. j2Ta)): Appar- 
ently the chosen spatial correlations in the disorder lead 
to a total breakdown of the TAI conductance plateau. 
We also observe an increased disorder-sensitivity as well 
as a narrowing of the conductance plateau correspond- 
ing to the original QSH phase in the clean limit. These 
results demonstrate that spatial correlations in the dis- 
order add an important new component to the physics 
of topological insulators. Especially in view of the en- 
visioned experiments that probe the physics of strongly 
disordered topological insulators, our results can appar- 
ently be expected to impose rather strict limits on the 
observability of the TAI phase. For such an experimen- 
tal realization of the TAI we can certainly conclude that 
it is equally important to be able to control the corre- 
lation length £ as it is to control the strength U of a 
disorder potential. 

In a next step we performed additional calculations to 
check the influence of spatial correlations also for larger 
samples than the ones considered above. This is particu- 
larly important as finite-size effects in small samples due 
to a coupling of counter-propagating edge stat es can co n- 
siderably distort the phase diagram of the TAl! 27 ' 28 ' 4 ^ To 
determine these phase boundaries in extended samples 
we performed a scaling analysis following previous work 
in this direction ! 19 ! 45 ! For this purpose a quadratic geom- 
etry of size L — L x — L y is rolled up to a cylindeJ 
using periodic boundary conditions in y-direction which 
eliminate the edge states in the sample. The disorder- 
averaged logarithmic conductance (In G) of the remain- 
ing bulk states is then calculated for three different sizes 
L\ < L2 < L3 of the quadratic system. The cross- 
ing points of (InG) calculated for samples of different 
sizes determine the disorder boundaries at which the 
bulk states make transitions between different charac- 
teristic transport regimes.^ We calculated these phase 
transition points for an uncorrelated and for correlated 
potentials with different correlation length £ in systems 
of three different sizes L\ — 500 nm, L2 = 700 nm and 
L3 = 1050 nm. Due to the high numerical effort involved, 
we restrict ourselves to a single energy of E F — 16 meV 
at which the TAI conductance plateau in the uncorre- 
lated case is wide and well established (see Fig. [2](a)). 

The results of our scaling analysis are shown in Fig. [3] 
Considering first the uncorrelated case studied already 
earlier (see black curves) we find that the lowest cross- 
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FIG. 2. (color online). Average conductance (G) as a func- 
tion of disorder strength U and Fermi energy Ef through a 
topological insulator (m = — 10 meV). The system length 
L x = 2000 nm and the height L y = 500 nm. The average 
is taken over (a) 200 and (b) 1000 disorder realizations. The 
green area shows where the average conductance (G) reaches 
a plateau at a single conductance quantum e 2 /h originat- 
ing from edge transport through the disordered system, (a) 
For uncorrelated disorder the TAI phase appears at energies 
Ef > \m\ and strong disorder. The blue lines show the phase 
boundaries of the TAI predicted by the effective medium the- 
ory which is based on the self-consistent Born approximation 
(SCBA) shown in Eq. jHJ. (b) Spatial correlations in the dis- 
order with correlation length £ = 23.45 nm destroy the TAI 
conductance plateau and narrow the plateau to an interval 
within the original bulk band gap. 



ing point occurs at a value of the disorder strength 
U « 70 meV which fits well with the onset of the afore- 
mentioned TAI conductance plateau in Fig. 2|a) (see also 
Ref. [45]). This onset coincides here with the opening of 
a bulk band gap which is reflected in the scaling plots of 
Fig. [3] through a dramatic reduction of the conductance 
by more than ten orders of magnitude. The second and 
third crossing points, in turn, can be associated with the 
breakdown of the TAI phase observed in Fig. [2ja) . The 
corresponding transition is, however, not induced by a 
band edge but rather by a mobility edge associated with 
those bulk states that fill the band gap when increas- 
ing the disorder strength beyond the first crossing point. 
These bulk states undergo a delocalization-localization 
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transition at the second and third crossing point which 
destroys the conductance plateau as soon as the delocal- 
ized bulk states start coupling the edge states at opposite 
edges in the sample. Note that for this to happen it is 
already sufficient for individual rather than for all bulk 
states to delocalize such that finite size effects d o play a 
role at this strong-disorder boundary of the TAl ! 45 * 46 J 

When extending the above scaling analysis now to the 
case of correlated disorder with successively increasing 
correlation length £ we find a behavior different from 
the uncorrelated case: already for the case of very short- 
range correlations with £ = 9.0 nm (red curves in Fig. [3]) 
significant differences appear. We still find three crossing 
points as before, but the conductance no longer displays 
the very strong suppression associated with a bulk band 
gap. Instead, we find that the delocalization-localization 
region of bulk states, which was previously associated 
with the strong-disorder boundary of the TAI, widens for 
increasing correlation length £. Correspondingly, in the 
disorder interval between the first and the second cross- 
ing points (which are now also much closer together) the 
conductance is much less suppressed than in the uncor- 
related case. This indicates that in the correlated case 
the bulk band gap disappeared and was filled with local- 
ized bulk states. We emphasize, however, that both a 
band gap as well as localized bulk states can give rise to 
a TAI, as was explicitly pointed out in a recent studyP^ 
In what was termed a TAI-I phase the coupling between 
edge states is prevented by a bulk band gap which elimi- 
nates all bulk states that would mediate such a coupling. 
A second TAI-II phase was characterized by a coexis- 
tence of localized bulk-states and extended edge-states. 
As long as the localization length of these bulk states re- 
mains smaller than the width of the sample, the coupling 
of edge states remains suppressed and the TAI persists. 
Following these arguments, the transition into the TAI 
phase can either occur at a band edge (for TAI-I) or at 
a mobility edge (for TAI-II) . From this we conclude that 
the original band edge at the weak-disorder boundary of 
the TAI-I phase in the uncorrelated case gets replaced by 
a mobility edge as the new weak-disorder boundary of a 
TAI-II phase in the correlated case. The corresponding 
suppression of the TAI-I phase already at a rather small 
correlation length of £ = 9.0 nm suggests that the TAI-II 
phase is more robust to spatial disorder correlation than 
the TAI-I phase. The TAI-II phase is, in turn, more 
sensitive to finite size effects due to individual localized 
bulk states which can couple counter-propagating edge 
states. Correspondingly we can understand the absence 
of a TAI conductance plateau in Fig. [2jb) and Fig. [1] in 
between the first and second crossing points of our scaling 
analysis in Fig. [3] as a finite-size effect which may disap- 
pear for much larger samples than studied here. Further 
explicit calculations will also be necessary to better un- 
derstand the infinite-size limit for TIs with long-range 
correlations in the disorder. Our own results for the case 
of £ = 35 nm (see the blue curves in Fig. [3]) show that 
the widening of the delocalization-localization transition 
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FIG. 3. (color online). Average logarithmic conductance 
(In G) as a function of disorder strength U for different sys- 
tem sizes Li = 500 nm, L2 = 700 nm and L3 = 1050 nm 
for uncorrelated disorder (black) and correlated disorder with 
correlation length £ = 9.0 nm (red) and £ = 35.0 nm (blue) at 
an energy of Ef = 16 meV (disorder average taken over 2000 
configurations). The transitions into and out of the TAI phase 
can be seen by the crossing points of curves with equal color. 
In the correlated case the TAI bulk band gap between the first 
and second crossing point disappears and the delocalization- 
localization region between second and third crossing point is 
broadened. 



continues for increasing correlation length £. However, 
since the transition region for £ = 35 nm is here already 
very wide, detailed statements on the phase boundaries 
in the infinite size limit are difficult to deduce from our 
finite-size calculations. 

While further calculations will thus be required to elu- 
cidate the infinite size limit of the TAI-II phase, we can 
further corroborate our arguments from above with re- 
spect to the suppression of the TAI-I phase due to cor- 
relations in the disorder. If these arguments are correct, 
the QSH phase (characterized by a negative topological 
mass and a chemical potential in the band gap, < |m|) 
should border, in the correlated disorder case, directly on 
the surviving TAI-II phase (characterized by a negative 
topological mass and a chemical potential outside of the 
band gap, > \m\). In the corresponding plot in Fig. 
[2Vb) we see that the QSH conductance plateau (exist- 
ing in the energy-range —10 meV < \l < 10 meV in the 
disorder-free sample) extends only to much smaller val- 
ues of disorder strength U than in the uncorrelated case 
(compare with Fig. [2]Ja)). Since for the above arguments 
the borders of this reduced QSH conductance plateau 
with the neighboring TAI-II phase are characterized by a 
band edge crossing, these borders should be describable 
in terms of a similar effective medium theory as has been 
developed for the uncorrelated caseP^ This theory maps 
the disordered system onto a disorder-free sample with a 
renormalized topological mass m = m + 5m and chemical 
potential p, = fi + 8fx. This renormalization was carried 
out in terms of the self-consistent Born approximation 
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(SCBA) using an integral equation for the self-energy E: 

U 2 ( a \ 2 [« 
E = — ( — ] lim / dk T dk v 

x (E F + iK-Ho(k)-E) 1 . (6) 

Whenever the renormalized chemical potential reaches 
the edge of the band gap (\p.\ = \fh\) the border of the 
QSH phase has been reached. Since this indicator, as 
calculated through the above SCBA, is independent of 
the system size, the effective medium theory offers an 
insightful and practical tool to determine the boundaries 
of the QSH phase in the infinite-size limit. 

In order to generalize the effective medium theory from 
above to the case of spatially correlated potentials we re- 
sort to recent theoretical work in which an extension of 
the coherent potential approximation to correlated disor- 
der was proposed.^ Following this line of work, one can 
conveniently include the disorder correlations through 
an additional term in Eq. (JgJ) which is given by the 
normalized Fourier transform of the disorder correlation 
function^ <5(ifc), 

U 2 . [> 
E = — lim / dk x dk y 

12 K-S-0 J_H 

C{%) x (e f + in - H (k) -e) 1 . (7) 

Since in the present case of Gaussian disorder correla- 
tions [see Eq. Q] the expression for C(k) is a Gaussian 
itself (centered around k = 0), the effect of the disor- 
der correlations is to smoothly cut off the above integral. 
The more long-range the correlations are (in real space), 
the sharper this cut-off is (in Fourier space) . 

To extract the corrections 5m and S/i from the self- 
energy in Eq. |7|), we proceed along the lines of the un- 
corrected case^ and decompose the self-energy into the 
Pauli matrices Oi 

E = E ao + EjeCt^ + T, y (Ty + T, z a z (8) 

with the help of 

5m — Re E z , 
Sfj, = -Re E . (9) 

With the values extracted for the renormalized topolog- 
ical mass m = m + 5m and chemical potential /2 = 
Ep + <5/i, we can now estimate the boundaries of the 
QSH phase in the case of correlated disorder by determin- 
ing the values of Ep and U for which the renormalized 
chemical potential /2 drops into the effective band gap at 
p. = ±m. We start by first testing our approach for the 
uncorrelated case, for which the weak disorder bound- 
ary of the TAI phase was estimated before. In this case 
the discrete Fourier transform of the correlation function 
C(k) is constant. By normalizing this function in fc-space 




FIG. 4. (color online). The average conductance (G) as a 
function of Fermi energy Ef and disorder strength U is shown 
for the correlated case with correlation length (a) £ = 9 nm 
and (b) £ = 23.45 nm. The system is L x = 2000 nm long 
and L y = 500 nm high while the grid spacing a = 5 nm. The 
blue curves demonstrate the validity of our adapted effective 
medium theory from Eq. Q. 



to the volume V = (2ir) 2 /a 2 of the Brillouin zone, we ex- 
actly reobtain the expression for the self energy in the 
uncorrelated case, Eq. If we determine with this ap- 
proach the phase boundaries of the QSH phase in the 
uncorrelated case we obtain the blue curves in Fig. [2](a) 
which fit nicely to the conductance plateau of the Q SH 
and TAI phases and to previous calculations! 16 1 17 1 4 -^ Ex- 
tending our calculations to the case of correlated disor- 
der, the borders which we calculate through Eqs. Q and 
^ (see blue curves in Fig. [4]) describe the boundaries 
of the QSH conductance plateau very well (without any 
adjustable parameters). The good agreement which we 
find for different correlation lengths £ (see Fig. [4] parts 
(a) and (b)) corroborates the validity of our approach. 
Note that, in contrast to the uncorrelated case, no con- 
ductance plateau is observed for > |m| (i.e., outside 
of the energy region where the QSH phase is present in 
the clean sample). An increasing disorder strength rather 
leads to a narrowing of the bulk band gap within which 
pure QSH edge transport can occur. This reduced band 
gap corresponds to positive corrections 5m and 5[i in the 
case of a correlated potential whereas in the uncorre- 
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lated case these corrections were shown to be negative.^ 
It is exactly these positive corrections which lead to the 
breakdown of the TAI-I phase for > \m\. 

IV. CONCLUSION 

In this work we investigate the effect of spatially corre- 
lated disorder on a two-dimensional topological insulator. 
We thereby extend previous studies in which only un- 
correlated disorder potentials were considered! 16 * 17 * 4 ^ * 4 ^ 
Our calculations show that a finite correlation of the dis- 
order potential enhances finite-size effects and may en- 
tirely suppress a regime of quantized conductance known 
as the TAI phase. Such a suppression of the TAI phase 
has also been found for off-diagonal bond disorder origi- 
nating from the deformation of the latticed To describe 
the shifted boundaries of quantized conductance theoret- 
ically, we perform a scaling analysis and adapt an existing 



effective medium theory to the case of spatially correlated 
potentials which yields quantitative agreement with our 
numerics. Our results suggest that for observing the TAI 
phase experimentally, it will be necessary to work with 
comparatively large samples (to suppress finite size ef- 
fects) and with very short ranged disorder potentials as 
any long-range correlations may strongly suppress this 
topologically non-trivial phase. We speculate that spa- 
tial correlations might also be an important impediment 
to eliminate the bulk conductance in three-dimensional 
topological insulators.^! This would certainly constitute 
an interesting topic for further investigations. 
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